Interacting electrons with spin in a one— dimensional dirty wire connected to leads 
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We investigate a one-dimensional wire of interacting elec- 
trons connected to one-dimensional noninteracting leads in 
the absence and in the presence of a backscattering poten- 
tial. The ballistic wire separates the charge and spin parts 
of an incident electron even in the noninteracting leads. The 
Fourier transform of nonlocal correlation functions are com- 
puted for T > w. In particular, this allows us to study the 
proximity effect, related to the Andreev reflection. A new 
type of proximity effect emerges when the wire has normally 
a tendency towards Wigner crystal formation. The latter is 
suppressed by the leads below a space-dependent crossover 
temperature; it gets dominated everywhere by the 2k f CDW 
at T < L5( K - 1 ) for short range interactions with parameter 
K < 1/3. The lowest-order renormalization equations of a 
weak backscattering potential are derived explicitly at finite 
temperature. A perturbative expression for the conductance 
in the presence of a potential with arbitrary spatial extension 
is given. It depends on the interactions, but is also affected by 
the noninteracting leads, especially for very repulsive interac- 
tions, K < 1/3. This leads to various regimes, depending on 
temperature and on K. For randomly distributed weak impu- 
rities, we compute the conductance fluctuations, equal to that 
of 1Z — g — 2e 2 /h. While the behavior of VarfJZ) depends on 
the interaction parameters, and is different for electrons with 
or without spin, and for K < 1/3 or K > 1/3, the ratio 
Var(lZ) /1Z 2 stays always of the same order: it is equal to 
Lt/L <C f in the high temperature limit, then saturates at 
1/2 in the low temperature limit, indicating that the relative 
fluctuations of JZ increase as one lowers the temperature. 

72.f0.-d, 73.40.Jn, 74.80.Fp 



I. INTRODUCTION 

One-dimensional quantum wires provide an interest- 
ing opportunity to study mesoscopic physics in an in- 
teracting system. On the one hand, it is well estab- 
lished theoretically that interactions give rise to unique 
electronic properties in one dimension, described by .-the 
so called Tomonaga-Luttinger liquid (TLL) model.EJ A 
detailed comprehension of the remarkable interplay-, be- 
tween interactions and disorder has been achieved.eTEl On 
the other hand, little attention has been paid to the role 
of the contacts which are known to influence strongly the 
transport properties of mesoscopic structures and quan- 
tum wires. In this respect, two simplified models were 



proposed recently. Either one cojonects a finite wire to the 
reservoirs by tunneling barriers^ leading to the suppres- 
sion of the ballistic conductance at low temperature due 
to the dramatic effect of the interactions. Alternatively, 
in the opposite situation, the wire is perfectly connected 
to one-dimensional leads, LIB yielding a perfect conduc- 
tance independent oruthe interactions of any range less 
than the wire length. lI The latter result contradicts the 
conductance reduction,bv the interactions predicted in a 
wire without contactsJjB and is in agreement witipjecent 
experiments on micron-length quantum wires .I11JO The 
perfect conductance can be explained through an extep* 
sion of Landauer's approach to the interacting wire.E3 
The reservoirs are taken into account by the flux they 
inject which acts as an initial condition for the equation 
of motion of the density.Q The incident flux is perfectly 
transmitted for any range of interactions. But reservoirs 
inject electrons which are not the proper modes of the 
wire. This leads us to introduce intermediate noninter- 
acting leads so that we can properly identify the injected 
and transmitted electronic flux. Under theses circum- 
stances an incident electron undergoes multiple internal 
reflections at the contacts due to change in interactions, 
leading to a perfect, transmission into a series of spatially 
separated charges .0 

While the ballistic conductance of a quantum wire can- 
not reveal the TLL character of the wire, the natural 
question one asks is if the other spectacular manifesta- 
tions of the TLL model, that are the signature of its 
non-Fermi-liquid behavior, can be still observed, ft is 
the purpose of this paper to study these features .—Some 
have already been addressed in the literature,E3llj but 
new interesting effects of the leads will be seen to ap- 
pear. We salLalso discuss recent experiments on quan- 
tum wires JL3U 

Beyond these points, in the present paper we develop 
a formal framework to deal with the properties of any 
inhomogeneous Tomonaga-Luttinger liquid (1TLL). An 
ITLL can occur in circumstances more general than con- 
tact effects. For example, a spatially varying effective 
interaction can be due to a varying width of the wire 
or to a nearby gate with a peculiar geometry. One can 
also think of another ideal system, the edge states in 
the fractional quantum Hall effect (FQHE).Lll The ITLL 
model might be relevant to describe transitions between 
edges at differpiOillmg, or an edge state connected to 
a Fermi liquid.E£rE3 This motivates us to investigate not 
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only electrons with spin but also spinless electrons with 
external interacting leads. 

Let us now summarize our main results. A typical fea- 
ture of a TLL is the separation of the charge and spin dy- 
namics. Imagine that one injects a spin polarized flux of 
electrons into one external lead, and detects transmitted 
spin and charge on the second external lead. Since the 
latter is noninteracting, one might suspect that charge 
and spin recombine. But we show that this is not the 
case: the charge and spin parts are separated even in the 
noninteracting leads. 

Such process has been first studied in refs. f7||l5| for 
spinless electrons. The propagation is defined in terms of 
a quasiparticle current (corresponding to Laughlin quasi- 
particle in the FQHE), and to a superposition of electron- 
hole excitations in a quantum wire. The basic building- 
block in the description of this effect is the scattering ma- 
trix of one quasiparticle at the contact between a TLL pf 
parameter K\ and a second one of parameter K 2 .DO tlil 
The reflection coefficient is 7 = (K\ — K2)/{Kx + K2), 
thus is negative when K x < K 2 -r- fpich exotic result in- 
dicates an analogy with AndreevErEEl reflection even for 
repulsive interactions. 

In connection with this Andreev reflection, we are here 
led to show, a second aspect, the manifestation of proxim- 
ity effectsEOE3 It is known that a typical feature of the 
TLL is the nonuniversal algebraic decay of correlations of 
different types, reinforcing the tendency towards a CDW 
or superconducting order which cannot be of long range 
due to the importance of quantum fluctuations in one di- 
mension. It is important to know the effect of the finite 
size and of the leads on the fluctuations. We will show 
here that the dominant tendency in the wire extends to- 
wards the external leads. 

A new aspect of the "mutual proximity effect" man- 
ifests itself when the wire has a tendency towards a 
Wigner crystal (Akp^CDW). This happens usually for 
Coulomb interactionsE.3 but also for very repulsive short 
range interactions. For electrons with spin with charge 
interaction parameter K < 1/3, the Wigner crystal dom- 
inates the Ihp CDW at any temperature in an infinite 
wire, while the inverse holds at K > 1/3. With external 
leads connected to a wire with parameter K < 1/3, we 
show the existence of a crossover temperature depending 
on position T c (x) below which the 2kp CDW dominates 
the Wigner crystal. T c (x) has a non-trivial powerlaw de- 
pendence on both the distance to the contacts and the 
wire length. 

The behavior of density correlations can be probed 
through coupling to a backscattering potential. The ef- 
fect of weak impurities on the conductance of the wire 
with leads was ipwid to be qualitatively similar to that 
without leadsJl3L3iiJ't3 and this seems in agreement with 
experiments.Ej But compared to the usual TLL, the con- 
ductance is affected in a nontrivial way by the contacts, 
especially when isolated barriers are considered. Ejll3 We 
show here that for electrons with spin this dependence 
becomes more crucial for very repulsive interactions be- 



cause, as mentioned above, the external leads obscure the 
tendency towards the formation of a Wigner crystal at 
low temperature. 

Apart from the specific model we consider, we improve 
the study of weak impurities. Usually, the jjenormaliza- 
tion equations in the presence of one barrierO are derived 
for both the thermal length, Lt — Vf/T, and the wire 
length, L, infinite: their finite values are accounted for 
semi-empirically by introducing inf(£, Lt) as a cutoff. 
Here, the renormalization equations are explicitly derived 
for finite lengths L and Lt- Besides, they are extended to 
a backscattering potential whose total spatial extension 
is less than Lt and L. 

Another important question concerns the fluctuations 
of the conductance in the presence of a random potential 
distribution. This question can arise in experiments even 
if the mean free path is much larger than the wire length 
due to the imperfection and roughness on the boundaries. 
Here we show that the conductance is self-averaging at 
high temper ature,pinore precisely Var(TZ)/TZ 2 is of the 
order L/Lt <C 1-lJ But this ratio saturates at 1/2 at 
temperatures low compared to T& = u/L. This has to 
be taken seriously into account when one tries to infer 
precise powerlaws from the experimental value of g. 

Let us now give the plan of our paper which does not 
follow the above sequence of results. Instead, we sepa- 
rate the paper in two parts: the first deals with spinless 
electrons which are also relevant for the FQHE, while we 
restore spin in the second part. For spinless electrons, 
we begin by describing the formalism and the inhomo- 
geneous Luttinger liquid model. We discuss the discon- 
tinuous interaction case which is suitable to understand 
qualitatively the results; it might be relevant for edge 
states in FQHE, but is has to be taken with some cau- 
tion for quantum wires. 

We then study the density-density correlation func- 
tions for two reasons. First, in order to show how the 
dominant tendency in the wire extends towards the ex- 
ternal leads, in analogy with proximity effects. Next, to 
investigate the role of backscattering. In particular, the 
nonlocal correlation functions at finite temperature and 
at low frequencies are computed. Any extended non- 
random weak potential is then considered, for which the 
leading renormalization equation is derived explicitly at 
finite temperature. The conductance is expressed per- 
turbatively in the backscattering potential. We discuss 
the random-impurity case where the most important new 
result concerns the conductance fluctuations. 

In the second part dealing with electrons with spin, 
most of the spinless treatment can be extended except 
for two important points. First, the transmission pro- 
cess is studied explicitly. Secondly, we show in which 
way leads affect the tendency towards the Wigner crys- 
tal that exists for very repulsive short range interactions. 
This induces different regimes for the correction to the 
conductance in the presence of a backscattering poten- 
tial. Finally, we discuss the experimental results obtained 
in quantum wires. 



2 



II. SPINLESS ELECTRONS 
A. Model 

In a strictly one-dimensional system with short range 
interactions between spinlcss electrons the low-energy 
properties can be fully parameterized by two parameters 
u and K. Without interactions u = vp and K = 1. The 
Hamiltonian can then be expressed as 



[1=1* 

2vr 



uKU 2 + — (d x $ 
K 



(2.1) 



where $ is related to the long-wavelength component 
of the electron density through p = — d x $/ir, and II is 
the canonical momentum conjugate to <E>, [<i>(x), n(y)] = 
iir5(x — y). Taking into account the discreteness of the 
electrons, Haldane derived the representation of the total 
electron density in terms of 



E 

m=— oo 



(2.2) 



where $(x) = $(x) + kpx. The coefficients c m obey 
c_ m = c m . For the noninteracting case only cq and 
c±\ are nonzero, however in the presence of interactions 
higher coefficients appear. The c m then can be calculated 
perturbatively, as done implicitly in ref. 

Consider now an interacting wire delimited by [—a, a] 
whose length will be denoted by 

L = 2a, 

connected perfectly to non-interacting lea ds. The global 
system is described by the Hamiltonian (2j|) with spa- 
tially varying parameters u(x), and K(x)i!\'Q We require 
u and K to be uniform on the external leads, i.e. for 
| a; | > L/2, taking values ul and Kl- Even though 
the most relevant situation of noninteracting leads cor- 
responds to u = vp and Kl = 1, it is interesting to 
keep -Kpfc-jfor other possible applications, for instance the 
FQHEo 

We note that the absence of translational invariance in 
the quantum wire gives rise to an inhomogeneous chem- 
ical potential p(x) that can be incorporated by a trans- 
lation in <FE3 Then one has to replace, in eq.(|2.2|), 



kpx — > 4>o(x) = —kpx + 



/dx fi — . 
u 



(2.3) 



Let us discuss the microscopic arguments for the ITLL 
model. If one assumes that the screening does not 
take place in the same way in the two-dimensional 
gas into which the wire opens and inside the wire we 
can argue that interactions are described by a func- 
tion U(x,y) which is not translationally invariant. But 
then the electronic momentum is not conserved any- 
more in the vicinity of the contacts. Indeed, this ap- 
pears when one bosonizes the interaction Hamiltonian 



/ J U(x,y)p(x)p(y), using p from eq.(2.2): the expo- 
nential terms in <E> cannot be ignored in general, and 
those terms violate momentum conservation. But if 
U(x,y) = f(x)f(y)h(x — y) where / varies slowly on 
scales \p. one can reduce the Hamiltonian to a quadratic 
form l 16 H 12 l Strictly speaking, this condition is required to 
find a perfect conductance, as discussed in ref. 12. If this 



condition is not met, and for not too strong interactions, 
the perfect conductance is recovered in the high tem- 
perature limit, but is reduced due to the abrupt change 
in interactions at low temperature. For instance, when 
u(x) = u and K{x) = K are constants for |x| < L/2, 
exponential terms in <E>(±a) cannot be ignored, and give 
rise to effective barriers determined by the jump of the 
interactions. Nevertheless, due to the symmetry of the 
structure, and since one needs kpL 1 to be able to use 
the TLL model, one can easily achieve resonances that 
suppress the role of these terms. 

For quantum wires opening adiabatically into a two- 
dimensional gas, the situation of abrupt variations is not 
realistic; it is adopted for mathematical convenience, and 
the results can be be trusted only far from the contacts 
compared to where the behavior should not depend 
on the profile of variations at the contacts. We shall 
however look at the behavior at the contact because it 
is very similar to that in thcj-strong tunneling limit of 
two coupled Luttinger liquids:E3 the effective .parameter 
at the contact, yielding the locajLconductanceQ and con- 
trolling the correlation functional is exactly the same as 
the parameter for the equivalent TLL liquid in ref. 

The situation for edge states in the fractional quantum 
Hall effect (FQHE) might be different. Non-quadratic 
terms arise from tunneling between edge states, which is 
more difficult to achieve due to their spatial separation. 
Recently, the model with discontinuous parameters was 
shown to be suited tj&.describe transitions between edges 
with different filling ,t3 opthe connection between an edge 
state and a Fermi liquid.E3 In any case, most of the results 
we derive here can be extended to any profile of variations 
of the parameters. 



B. Correlation functions 

It is well known that symmetry breaking phase transi- 
tions cannot occur in one dimension. But the interactions 
in an infinite wire enhance charge density or supercon- 
ducting fluctuations depending on whether they are re- 
pulsive or attractive. The natural question is to know 
whether such tendencies persist in the presence of leads. 

In the present subsection, we would like to compute 
the correlation functions not only to ansAwetptiiis question, 
but also to explain the proximity effect ,Ercjll3 and to give 
the necessary tools for studying the role of backscatter- 
ing in the next section, closely related to the tendency 
towards formation of a charge density wave. 
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1. Correlation functions at finite temperature 

Let us wri te the density-density correlation function, 
using eq.(2.2): 



(T rP (x,T)p(y,0)) = -d x d y 



U(x,y, 



y -^X m (x,y,T)e 2lm[Mx) - My ^ 
2m 2 



(2.4) 



where U is gi'ven in terms of the fundamental bosonic 
Green functioned 



G(x,y,r) = (T T $(x, r)$(y, 0)) 



(2.5) 



as 



U{x, y, t) = G(x, x, to) + G(y, y, r ) - 2G(x, y, r). (2.6) 

To is a cutoff different but in general of the order of 1/Ep. 
We impose on the imaginary time r > To, so that G can 
be computed at finite temperature, without dependence 
on r .EaThe consequence of this cutoff procedure on an- 
alytic co ntin uation is elucidate d in appendix |^. 
In eq.([2.4|) (p is given by eq.(2.S), and 



im$(x,T) — im<&(y,0) 



X m (x,y,T) = (T T e % '--^'-'e 

,-2m 2 U(x,y,T) 



(2.7) 



Note that the first term (in U) in eq.( |2.4| ) can be obtained 
from the limit m — ► of X m . 

Instead of writing out explicitly X m the imaginary part 
of its Fourier transform in the low frequency limit to <C T 
will be of more use: 

lmX m (x,y,(j —r=^Xm{x,y), (2.8) 

F 

where Xm(x, y) is a dimensionless function. The function 
Xm can be obtained from xi by multiplying K(x) by m 2 , 
thus we will often study Xu denoted x f°r simplicity. 
The time dependenc e of U being due to the nonlocal 
part G(x,y,r) in eq.(2.6), we can write (for details see 
appendix H) 



X (x,y) = -^C(x,y)e- 2G ^ ) e -2G(y, y ,r )^ 



T 



where 



T = Trn 



(2.9) 



(2.10) 



and C is a dimensionless time integral, eq.(B9). We can 
find the properties of C(x,y) for general smooth varia- 
tions of K and u, and we have computed it explicitly for 
any x, y in the discontinuous interaction case and at finite 
temperature. By deforming the contour of integration in 



the complex plane, the computation is made much eas- 
ier, especially in the low temperature limit. In the high 
temperature limit, more steps are needed. The results 
are given explicitly in appendix [b| 

In an infinite wire with uniform parameter K, x{ x iy) 
is a function of x — y, and we denote it by x \ x ~ 
y). It decreases exponentially at separations larger than 
L T = ttw/T, x (K) \r) = {r/L T )e- Kr / LT , while its local 
value has the typical powerlaw in temperature X^ K \ r = 
0) ~ 7 i 2(A"-i) faofc diverges in the zero temperature limit 
when K < 1, indicating an enhancement of the density 
fluctuations compared to the noninteracting case. 

For interaction parameters u, K constant in the bulk 
of the wire and reaching their asymptotic values ul,Kl 
on the external leads, four general properties of x can be 
shown. The first and second one hold at high tempera- 
ture, where comparison can be made with X^ K ^ ■ 



1. At T > T L , we have 



a- M,a - \y\ > L T x{x,y) ^ X K x -y) 

(2.11) 

where x^ K \ x ~ y) 1S the value of x m an infinite 
wire with parameter K. Thus we recover the trans- 
lationally invariant behavior for any x, y far from 
the contacts compared to Lt- 

2. Again, at T S> Tl, but for any location of x, y, we 
have 



x( x ,y) < x {K) ( x -y) 



(2.12) 



for Kl > K, while the inverse inequality holds for 
K L <K. 

3. At any temperature, and for small separations com- 
pared to Lt, X factorizes: 



\x-y\ <L T 



X{x,y) ~ y/x( x ,x)Vx(y,y)- 
(2.13) 



Thi s is due to the fact that C(x,y) ~ C(x,x) 
[eq.( 2.9)1, U P to corrections of order (\x — ?/|/Lt) 2 - 



4. The local value of x is almost constant on segments 
much less than the wire length and the thermal 
length. If we denote 



then 

\x-y\<L 



Lmin = min(L, Lt), (2-14) 

=^ x( x , x ) - x(y, y) ~ x( x , y)- 



(2.15) 



These equalities hold up to corrections of the order 
of (|a; - y\/L min ) 2 . 
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Let us now specialize to the case of discontinuous pa- 
rameters. First we introduce the quasiparticle density 
j r = j + rup for r — ± which are the right- and left- 
going proper modes in a uniform TLL. They correspond 
to Laughlin quasiparticle currents in edge states of the 
FQHE. As we said before, the basic quant ity n £ojj— un- 
derstanding the physics is the scattering matrixQliXL3 of 
quasiparticles at the contact of two perfectly connected 
TLL with parameters K and Kl- It acts in the space 
(j+,3-) and is given by 



S 



—7 I + 7 
1-7 7 



where 7 is the reflection coefficient 
Kl - K 



7 



K L + K' 



(2.16) 



(2.17) 



The matrix S relates the outgoing flux to the injected 
flux, defined not in term of wavefunction amplitudes as in 
usual scattering approaches, but directly in terms of the 
quasiparticle current j r . The locaL-jeflective parameter 
(yielding a local Kubo conductanceBEa) iSj-given by K\ 
multiplied by the transmission coefficient ,0o 



K a = K L (\-i)=K(\+i). 



(2.18) 



Let us first discuss the high temperature limit. It is 
worth noti ng tha t a lower (upper) bound to x(x, y), com- 
pleting eq.( 2.12 ), is g iven by x^ K { x ~ y) f° r Kl > K 
(Kl < K). In cq.(|2.9D, the nonlocal part C(x, y) is given 



by eq.(B14) in the l imit T ^s> Tl, while it simplifies to 
a constant |eq.( |B20| )l in th e limit T Tl, in accor- 
dance with prop erty ( [2.13 ). The Green function G is 
given by cq. flBlO] ). Instead of wri ting explicitly the fac- 
tor e~ 2G ( x ' x ' r °) that enters in eq.( |2.9| ), we use the local 
value of x- This is equivalent as far as the dominant be- 
havior is concerned because C(x,x) is a slowly varying 
function of x (cf . appendix Ib|) . Thus we drop it from 



-2G(x,x,t ) 



(2.19) 



In the high temperature limit (T > Tl), the multiple re- 
flections caused by the change in interactions don't affect 
the correlation functions due to the lack of thermal co- 
herence along the wire. Only one reflection is felt within 
a thermal length Lt — u/T from the contacts. We now 
define 



a — \x\ 

UTq 



(2.20) 



the time it takes for a quasiparticle emanating at x to 
reach the closest contact, measured in units of To. For 
t x > 1 [eq.(|2.2CD], we find 



y/x{x,x) ~T K -\tanhTt x ) 



\K a —K 



(2.21) 



Again, for points far from the contacts compared to Lt, 
this expression coincides with that obtained in an infinite 



i2(K-l) 



At the contacts, 



TLL with parameter K, x ~ T_ 
we find x(a,a) ~~ jV{K a -i) j rf eq (pp^j ( this amount s 
to setting t x = 1 in eq. ( 2. 21 ), using T< 1). 

In the limit of temperatures low compared to Tl, the 
multiple reflections affect the correlation functions, thus 
the external leads with parameter Kl determine the tem- 
perature dependence. But there is a nontrivial depen- 
dence on the wire length as well as on the distance of the 
barrier from the contacts: 



VxM-Z^-'T^H 



K„—K 



(2.22) 



This gives the behavior for t x 3> 1. x( a i a ) can be ob- 
tained by setting t x = 1 in this expression. 

Remember that one obtains Xm(x,x) from x( x i x ) ( m 
both temperature limits) by mu ltiplying both K and Kl 
by m 2 , which leaves 7 [eq.( 2.17 )] unchanged. Clearly, the 
2mkF CDW is less important than the 2Uf one. Fur- 
ther, if electrons in the external leads are non interact- 
ing, Kl — 1, the 2mfci? CDW are suppressed in the zero- 



temperature limit because x-n 
[eq.(H)]. 



^2(m 2 K L -l) 



= 7 l2 ( r ' 



2. Proximity effect 

Let us first recall an exotic phenomena emanatuigJxDm 
this model, the analogy with Andreev reflectionEo'EEl 

a. Andreev reflection The reflection coefficient of a 
quasiparticle from a TLL with parameter K\ incident on 
another one with parameter K2 is given by 7, eq.(2.17), 
and is negative for K\ < K2, i.e. a partial quasi- hole is 
reflected back. This is similar to Andreev reflection of 
an incident electron on a normal-superconductor inter- 
face. The similarity is the closest for K\ = 1 and K2 > 1, 
because the wire has a tendency towards superconductiv- 
ity, but it holds more generally even for repulsive inter- 
actions, for instance for a quasiparticle in the wire with 
K < 1 incident on the noninteracting lead. 

It is worth noting that in the case of two semi-infinite 
wires with K\ and K2 as parameters, and K\ < K2, the 
local Kubo conductance at the interface is given by the 
transmission coefficient from 1 to 2 multiplied hy. Kx, i-e. 
by g a = K a . This verifies K x < g a < 2KiEEM3 This is 
very similar to the inequalities for an N-S interface, if one 
interprets K\ as the conductance (Kubo) of the normal 
side, g N < g N ^ s < 2g N SB 

Finally, let us comment on a recent paperEJ based on 
the model in ref. |2^ of two TLLs connected by a tun- 
neling term. The bosonized theory was associated with 
conformal field theories to describe quasiparticles in term 
of soliton states. Interestingly, in the strong tunnel- 
ing limit, the scattering matrix for these states at the 
interface turns out to be the same as eq.(2.16). This 
gives a stronger foundation of the quasiparticle scatter- 
ing scheme we introduced previously.Q 
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b. Proximity effect The reflection at the contact gives 
rise to a proximitg|-e|?efjt-|for both superconducting or 
CDW correlations.0H't3'Ell This can be seen computing 
the correlation functions on the external leads. In par- 
ticular, when T Tl, there is no coherence between the 
cndpoints of the wire, and it is sufficient to consider half 
the system. The wire and one lead have a sy mme tric 
role; it is sufficient to permute K and Kl in eq.(2.21) in 
order to find x for |a;| > a, and to take t x = {\x\ — a) Jul. 
By the way, this situation is relevant for the FQHE, if 
two edge states with differe nt fil ling are connected.!!^ 

At T < Tl, instead of eq.( 2.22 ) inside the wire, outside 
the wire one has 



^/^x)-T KL -\T L t x ) 



K a -K L 



up to t x TL ~ 1. Beyond a distance of the order of L in 
one external lead, we recover the simple law T Kh . 

For clarity, we discuss the consequence for nonintcract- 
ing leads, Kl — 1. Then if K < 1, the density-density 
correlations in the bulk are similar to those in an infinite 
wire, but are reduced when one goes to the contacts be- 
cause t]. K decreases for 7 > 0. Besides, this enhancement 
extends irir-the, external leads up to a distance L min /K 
[eq.(2.14)]£3'E3 For K > 1, the pairing correlation func- 
tion, which can be obtained from \ if one replaces K 
by 1/K and Kl by 1/Kl, is enhanced up to a distance 
KL m in- This increases with K, i.e. where interactions 
are mpxe attractive. This is reminiscent of the proximity 
effect.B 

In general, we can show that the proximity effect ex- 
tends up to the minimum length scale at hand. We must 
note that this holds even for a smooth profile of u and 
K, as for instance the properties of x exposed before are 
general. 



The forward scattering term — d x <f>V(x) is absorbed by a 
translation of $ by 



dx'—Vix 1 ) 
u 



(2.24) 



which has to be included in the scalar function 



eq. flgj) 



!>0. 



Note that the impurity Hamiltonian (2.23) is similar 
to that considered usually if we replace our V'(x) by 
—2imV(x) I (c m UTo). For instance, if V has the form of a 
kink, it corresponds to what is commonly treated as one 
barrier. 

In order to derive the renormalization equations, the 
exact partition function Z at finite temperature is ex- 
panded in terms of V: 



C. Backscattering by non— random impurities 

We study now the role of impurities in the wire. The 
conductance with impurities was shown to depend on 
the irtter|actkms and is generally affected by the external 
leads The main goal of this part is to give a gen- 

eral scheme to treat the effect of a weak backscattering 
potential V(x) in a non-translational invariant system 
and at finite temperature. 



1. Renormalization Equations 

The coupling of the conduction electrons to impurities 
wit h po tential V(x) is J dxV(x)p(x), where p is given by 



eq.(2.2). Itris then convenient to do an integration by 



parts giving! 



S2 



m— — oo,m^0 

(2.23) 



G 



(-1)" 



2J rmmjUij 
ft 



(2.25) 



where [Zy = [/ (xi,Xj,Ti — tj) is given by eqs.(2.5, 2.6), and the integration runs over Xi and T\. Z describes a neutral 



gas of integer charges restricted to a cylinder whose perimeter is (3, and whose height is deter min ed by the spatial 
extension of V, reducing to a circle when V' is local. The charges interact via U(x,u,t), eq.( |2.6| ). This is not the 
Coulom b int eraction but an infinite scries of logarithmic terms related to the transmission process discussed earlier 



(see eq.(B10)).E3 The renormalization procedure is implemented by increasing the cutoff to To(l) = tqc 1 , where tq is 



the bare cutoff 1/A, and modifying the para meters in order to keep Z invariant. It is immediate to derive the leading- 
order equation for V. According to eq.(2J3), the cutoff appears in the local part of U. The change of U(x,y,T) due 
to a change of the cutoff by cLtq is 



dr (l), 



The exponential of such a term for a pair Xi,Xj factorizes and, once the global neutrality is used, can be absorbed 
separately into V'(xi) and V'(xj). This gives the leading-order flow equations for V(x) explicitly at finite temperature 



dl 



V^(x,l) [\ + 2m 2< ^{x,x,T Q e l ) 



(2.26) 



Note that in the limit of zero temperature, an infinite uniform wire with parameter K, and a purely local V, t 
Fourier transform of this equation yields the well-known flow equation for V(2mkp) = V m , dV m /dl = V m (l — m 2 K)j 
provided kp is independent on the cutoff as we required. 
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In general, one expects V(x) to renormalize the inter- 
actions. In the extreme case of a local barrier, it was 
shown by integrating out degrees of freedom away from, 
the barrier that the interactions are not renormalized .El 
We can both recover and generalize this result in a dif- 
ferent way, by using the expansion above at finite tem- 
perature and in the finite wire: new interaction terms W 
between the charges are generated by the renormaliza- 
tion, but they decay faster than U at long time scales. 
For instance, in the zero-temperature limit, and for a bar- 
rier in the center of the wire, U = K log r is corrected by 
W(t,T = 0) = K' (log t)/t, where K' is renormalized by 
V. At finite temperature, we have W-(t, T) < W(t, 0), 
i.e. W is still decaying faster than [/.Il3 

Let us now discuss a more extended potential than 
one barr ier, but whose total extension obeys d <C L mm 
|eq.( 2.14 )]. It turns out that when one goes to a cut- 
off such that utq{1) > d, the partition function ( |2.25| ) 
becomes identical to th at of a local barrier.EEl This is 
related to the property ( 2.1 5| ) . Then it is tempting to 
assume that the effect of such a potential would be sim- 
ilar to that of a local barrier at length scales larger than 
the potential extension d. But one has to check that 
the renormalization from the bare cutoff up to d does 
not modify the local interaction parameter, i.e. K(x) at 
points where V'(x) ^ 0. This would induce non transla- 
tional invariant effective interactions even if the bare K 
is uniform. We can for instance show that tfio [eq^iO)] is 



V m = / dxV (x)e 2mMx \ 



(2.29) 



renormalized by a complicated complex function, whose 
equation is rather lengthy to write. Thus we think this 
point needs a more thorough study. 

Finally, it is worth writing the next-order corrections 
to eq.(2.26) in the case where the partition function be- 
comes equivalent to that of a barrier and therefore the 
integrations over Xi can be done explicitly in eq.(2.25). 
For given m, one has to add to this equation 



E 



(2.27) 



up to nonunivcrsal prefactors depending on toi,to 2 Ji1 

If we ignore renormalizations of the interaction, as is 
surely correct for a local barrier or for a weak enough 
extended potential, the equation (3.1C) can be integrated 
straightforwardly. In the absence of an external energy 
scale, the unique limitation on increasing the cutoff comes 
naturally from the fact that one has to put at least two 
charges on the same cylinder of radius (3, thus we stop at 
r = (3/2. Then the renormalized potential is simply 



(2.28) 



where xi x > x )> eq.( 2.19 ), has been studied in the previous 
subsection. 

One can also consider the Fourier transform V(k) of 
the potential. If th e extension of V is less than L m ; n , 
we can use eq.( 2.15| ) to factor out \ from the integral, 
thus allowing us to recover the dominant term of the 
renormalized Fourier component of V: 



by simple multi plica tion by ^Jxm{x,x). Note that 
whenever </>o [eq.( |2.3| )] can be replaced by kpx, V m — 
—2imkpV {^Imkp) . However the nonlocality of Xm has to 
be taken into account when one looks to the next leading 
term.Lj 

The derivation above is valid for any profile of parame- 
ters. In case K varies from its external value Kl towards 
a plateaus at K inside the wire, we can show that the de- 
pendence of x(cc, x) on T is monotonic. Its monotony is 
determined for any x by the sign of 1 — K at T > Tl, thus 
is similar to that in an infinite wire. The renormalized 
V , Eq.(2.28) obtained at a given temperature goes up 



(down) when temperature is lowered for K < 1 (K > 1). 
At T < Tl, the T dependence is controlled by the sign of 
1 — Kl, in particular V saturates for Kl = 1 at T < Tl- 
In the case of discontin uous param eters , the explicit 
value of x is given by eqs.(2.21) and ( [2.22 ) for the high 
and low temperature regime. The renormalized V, de- 
pending on the location of the barriers, is shown in fig.|l| 
for the case Kl = 1, K < 1. We point out that these 
curves are not inferred from cutting the scaling at T or 
at Tl, but result from the renormalization at finite tem- 
perature and for the finite wire. Note that the sat urati on 
at T < T L occurs only for K L = 1 and for m — 1 Jl3l2J 



2. Conductance with non-random impurities 

Let us first recall that the conductance of a wire with- 
out impurities g = e 2 /h is independent of interactions if 
the external leads are noninteracting. The conductance 
can be related to the transmission, which turns out to 
be perfect, thus-generalizing Landauer's approach to an. 
interacting wireLrE3 One can also use the Kubo formula,S 
provided the external leads are noninteracting, and one 
is limited to the linear and stationary regime. We have 
to stress an important point: the external field has to 
be used when interactions are included exactly in the 
Hamiltonian, and not the internal electric field modified 
by the interactions. This is clear for instance through the 
equation of motion method we have developed in ref. 
to study transport, for which the external field forms a 
source term. 

The same conditions to justify the Kubo formula hold 
with impurities, if one assumes that the reservoirs im- 
pose an external bias AV^ (that can be maintained in the 
noninteracting leads) independently of the current. Then 
the stationary current is given by 



j(x) = lim / cr(x,y,uj)E(y,uj) 
= a(x,y)AV 



(2.30) 



where a(x,y,Lo) is the nonloca l co nductivity, related to 
the Fourier transform of G [eq. (|2.5| )1 computed now with 
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-1.5ln(T L ) 




4 6 

-ln(T) 



FIG. 1. The three lower curves show the renormalized bar- 
rier strength, whose square yields the dimensionless reduction 
to the conductance TZ, scaled by its bare value. The dotted, 
dashed and dot-dashed curves correspond respectively to a 
barrier at the center of the wire, at the contact and at an inter- 
mediate point x. T is in units of 1/to. We chose to = Tl/148, 
thus — ln(TtTo) = 5, and K p — 0.5, K a — 1. The same curves 
hold for spinless electrons if one replaces K p by K, but one 
has to multiply the loglZ by 2. There is a crossover at Tl 
from a powerlaw controlled by K p (K ap ~ 2K P /(1 + K p )) at 
the center (at the contact) to a plateau. For a barrier at x, 
there are two crossovers, one at T x = u p /t x from T Kp ~ 1 to 
T Kap ~ 1 , followed by a saturation occurring at Tl- The con- 
tinuous upper curve corresponds to l e (TV) /(uto) for an ex- 
tended Gaussian distribution, where the powerlaw at T < Tl 
is governed by K p since the interactions are repulsive. The 
2&F backscattering dominates the AUf one for any T since 
K p = 0.5 > 1/3. 



the total Hamiltonian H + Hi mp , and continued to real 



frequency, 



v(x,y,u) = &{x,y,u) 

IT 



(2.31) 



The uniformity of the current as well as time reversal 
symmetry requiije.-the zero frequency limit of o~(x,y,uj) 
to be uniformity Thus 



9 = -r<r(x,y)- 



(2.32) 



If the external leads are interacting, one cannot impose 
an external value of the potential, rathcti-thc potential is 
now renormalized by the interactionsJaO Nevertheless, 
the spatial separation of the two edges in a Hall bar can 
lead to different couplings with reservoirs as discussed 
in [H],^0], and it is also possible to measure locally the 
potential on each edge. Thus it is still of interest to let 
the external parameter Kl arbitrary. 



In appendix we derive a novel exact Dyson equation 
for the conductivity in a non-translationally-invariant 
system. This is very suitable ±q write the perturbative 
expression for the conductance!!!] 



Qirnp 

where TZ is given by 
TZ = 



Kl 
2tt 



rn—l 



(2.33) 



(2.34) 



and TZ m is the contribution from backscattering of m 
electrons, 



n,. 



Xm(x,y) 



V'{x)V'{y) 

{uk F y 



cos2m [(j>o{x) - (po(y)] , ( 2 -35) 



where \m is defined by Eqs. (2.6,2.5jEq,2.7). Remember 
that the forward scattering contribution is included in (f>a 
[cq.( 2.24 )]. We have already obtained general properties 
of %(x, y). If parameters change abruptly, one can inject 
the explicit form (appendix pf) to express the conduc- 
tance at any temperature, for sufficient weak potential 
with any extension. 

If the potential extension is much less than L m .;„ 
[eq.(2.14)], we can use cq.( 2.15|) t o show easily that the 
dominant term of equation (2.35) reduces to the square 
of the renormalized potential 



Xm{x,x)\V m 



(2.36) 



where V m is given by eq.( 2.29 ). As long as Vi ^ (i.e. 
out of resonance), the dominant term in 1Z is given by 
TZ±. On resonance, the dominant term comes still from 
m = 1 or from m — 2;E3 this depends on the values 
of K and Kl, as well as the extension and strength of 
the potential. It appears that at low temperature, be- 
cause the 4kp CDW correlation function acquires a factor 
T 2(4K L -i) [ eq- (|2,22l) with K and K L multiplied by 4], one 
must have Kl > 1/4 for TZ to vanish at zero temperature 
(for a short enough wire). Remarkably, for Kl = 1, the 
low temperature conductance is still e 2 /h on resonance, 
even for K < 1/4. We will not discuss the resonance 
in further detail here, even though one has a variety of 
beha viors, an d we refer to the explicit expression of TZ, 
Eqs. (2.34 2.35] ), and to the computation of appendix^. 

Out of resonance, and for discontinuous parameters, 
Kl = 1, 72. is given by the same curves as those yielding 
the renormalized V in fig.|l|. The wire has to be short 
enough so that TZ stays weak enough when reaching the 
plateaus below Tl, otherwise the perturbative computa- 
tion breaks down at a temperature above Tl- We refer 
to ref. Eq for the strong barrier case. 







D. Random impurities: conductance fluctuations 

Now we consider random impurities distributed all over 
the wire. For simplicity, we limit ourselves to the case of 
a Gaussian distribution, (V(x)) = 0, and 



(V(x)V(y))=DS(x-y). 



(2.37) 



By averaging eq.( 2.35| ) over disorder, one obtains the 
leading term of TZ [eq.([T3J)] coming from the m = 1 
contribution 



(u) = k l r 



dx 



(2.38) 



where l/l e (x) — <jD\<j}' (x)] 2 / (uk p) 2 . The forward scat- 
tering part in O is KV(x)/u [cq.( [2.24 )], and thus can be 
dropped when averaging. Thus the dependence of l e on x 
is due to the inhomogeneity of the interactions. It could 
also be induced by a space-dependent disorder strength 
D(x). Nevertheless, we will assume 4>'q{x) ~ kp for sim- 
plicity, in which case l e is uniform and is equal to D/u 2 . 
This holds for instance in the bulk of the wire, and thus 
this is a plausible approximation whenever the impurities 
in the bulk dominate. 



Note that we have obtained eq. (2.38) by p erforming 
first a double integration by parts in eq.(2.35) and then 
retaining only the term with no derivative of x- The ex- 
plicit integration of the function \ can be done but i s te- . 
dious. However, the results can be understood easily£§t£l 
There are two main contributions: one from the bulk 
that behaves, for T > T L , as LT 2 ^ K ^ X \ the other from 
the contacts behaving as T 2 ( Ka ~ l \ p or y < j i Lj mJ# 
nonuniversal constants, these contributions becomdij'u 



(TZ) 



le 



I rp Z 



T 



2(K-K L )-1 



T, 



2{K a -K L ) 



(2.39) 



where K a is given by eq. ( 2.1§| ) . 

Let us discuss the case of noninteracting leads, Kl = 1. 
In this case, the contact contribution dominates for very 
attractive interactions, K > (3 + vT7)/2, and for not too 
high temperature s (or not too long wires). 

Inspecting eq.(2.39), we see that (TZ) saturates at 
T < T L at a value < 1 for K > 3/2, but that it in- 
creases with wire length for K < 3/2. In the latter case, 
the perturbative computation is valid only for L < L; oc 
where 



L 



loc 



UT 



(2.40) 



This coincides with the localization length inferred by 
scaling arguments in ||,[| 

Let us now compute the variance of the conductance 
Var(g) = (g 2 ) — (g) 2 . Clearly, it is also equal to the vari- 
anc e of TZ which can be expressed perturbatively, using 
eq.( 2.35| ), as 



/a po 
/ 
-a J —i 



dxdy 



cos 2 2 [<fo (x)~Mv)]x 2 (x, y). 

(2.41) 



Here we have used eq.(2.37) to perform the average of 1Z 2 , 
have dropped the product of two correlation functions for 
different pairs of inte gers, and neglected the role of the 
derivative in eq.( [2.23| ). 

It is interesting to consider first the low temperature 
limit T <Tl, where a general inequality, valid for arbi- 
trary profi le K (x), can be shown. Using the factorization 
property (|2.13D which holds now f or any x, y in the wire 
because L <c Lt, we can write eq.(2.41) as 



Var(K) = E ^ / 



r=±l 



ra dx 



X(x,x) 



(2.42) 



The term for r = 1 on the right hand side is less th an th at 
for r = —1, this latter is nothing but (1Z) 2 |eq.( 2.38 )1- 
Thus we have the interesting inequality 



\{K) 2 <Var{lZ)<{lZ) 2 . 



In particular, this shows that 



Var{K) _ 1 
(TZ) 2 ~ 2" 



(2.43) 



(2.44) 



If 4>q (x) is simply kpx, this is a very good e stim ate of 
Var(7Z) because the integral for r = 1 in cq. ( [2. 42 ) con- 
tains a rapidly oscillating function, e 2ltf-1 , while x varies 
slowly compared to n/kp. Besides, Lkp S> 1 in order 
to validate the bosonization, so that the integral corre- 
sponding t o r = 1 in eq.( |2.42 ) can be neglected. This 
justifies eq.( 2.44| ), showing that the fluctuations of TZ are 
of the same order as its average value. In the special case 
of short range interactions with discontinuous parame- 
ters, Var(lZ) is obtained as the square of eq.(2.39). 



We now consider the high temperature limit, T > Tl- 
To simplify the discussion, we restrict ourselves to the 
case K < 1, which is probably the case in a quantum 
wire. Then the bulk contribution dominates over the 
con tact effects, so that we can restrict the integral in 
eq.( |2.4l| ) to -a/2, a/2, and therefore t x T > 1 Ov er this 
segment, x depends only on x — y ( see cq.(2.11)), and 
we can write, after a change of variables x — y = vLt, 



Var {n) ~ ^T^ K - l h{2k F L T ) 



where 



no 



(l + cos2£u) [C K (v)] 2 dv, 



(2.45) 



(2.46) 



and Ck is given by eq.(B15). Since L» hp, we have 
extended the integral approximatively to infinity. But 
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since we are restricted to temperatures much less than 
the Fermi energy, we have £ = kpLx ^S> 1, thus the 
oscillatory part of I can be neglected and I tends to a 
constant. The ratio 



Var{K) Lt 
(ll) 2 ~ ~L 



(2.47) 



is now much less than unity, contrary to the low tempera- 
ture limit, thus the relative fluctuations of 1Z are small. It 



is worth noting that one can extrapolate the ratio ( 2.47 ) 
obtained for L ^> Lt to L ~ Lt, to recover the same 
order of magnitude as eq. ( flj) at which Var(K)/ (II) 2 
saturates at low temperature, i.e. for Lt > L. 



III. ELECTRONS WITH SPIN 

Let us take into account the spin of the electrons, which 
is necessary when dealing with quantum wires, and to 
confront the theory with experiments. This case might 
eventually also be relevant for two-channel edge states in 
FQHE. 



A. Bosonization 



The inhomogeneous TLL model is now characterized 
by x dependent functions u u (x) and K v (x). For sim- 
plicity, we take noninteracting leads, otherwise the pa- 
rameters would be too numerous; indeed this is the most 
relevant situation for quantum wires. Nevertheless, the 
treatment of g± is not trivial for non-translationally in- 
variant interactions because g± now is space-dependent. 
We hawe derived the renormalization equations for this 



case ,E-j nevertheless their integration is quite involved. 
We can however draw qualitative but not firm conclu- 
sions, restricting ourselves to spin isotropic interactions. 

The most difficult case corresponds to attractive inter- 
actions. As in the uniform case, we expect a spin gap to 
develop, but it is not clear how it extends spatially and 
how is it affected by the leads in the low temperature 
regime. The way K a would vary from 1 on the external 
leads to K* — inside the wire is not clear neither; this 
is an interesting problem to solve. 

For repulsive interactions, the inhomogeneous K a is 
expected to renormalize towards a value respecting the 
SU(2) symmetry K* = 1. Loosely speaking, this is eas- 
ier to study because the leads have K a — 1, so that 
they don't prevent this scaling and rather favor it. This 
seems the most relevant situation for quantum wires, 
even though we can sometimes allow for K a to take any 
value. 



The typical feature for interacting electrons with spin 
is the separation of the charge and spin degrees of free- 
dom at low energy. The density for each spin component 
p s (s = ± denoting the spin up or down) has a long- 
wavelength part that is related to a boson field $ s by 
p s = —dx^s/n. The total density of electrons with spin 
s can be written as eq.(2.2) where $ has now the sub- 
script s. The charge and spin fields are defined by 



The Hamiltonian describing the low-energy properties 
can be decoupled in a charge and spin parts, H = 
H p + H a , where 



H„= — 



dx 
2^ 



(3.1) 



for v = p,a. The boson fields are related to the charge 
and spin density (p and a) through ^/2d x ^ p (x)/7r = p 
and \/2d x $> (T /'K = a. 1I„ is the momentum density conju- 
gate to <1>i/. In the absence of interactions, K u = l,u v — 
v F . 

Recall that an additional term has normally to be 
added to H a [eq.(3.1)] corresponding to the backscatter- 
ing of electrons of opposite spin, g± J dx cos V8$ a . For 
spin-invariant interactions, g± renormalizes to zero if it 
is initially positive, while a spin gap is opened if it is 
negative. Then K a scales respectively to 1 or zero which 
are the only values consistent with SU(2) symmetry. 



B. The transmission process 

The first issue we consider is the transmission process 
of an incident electron. The case of electrons without 
spin was discussed in refs. 0,|l3|, and this can be extended 
easily to electrons injected from reservoirs. 

Let us imagine that we inject an electron of definite 
spin in the left lead. This corresponds to create a kink in 
both $ p and <& a whose time evolution has to be solved 
given the initial conditions. The equations of motion for 
these fields are decoupled, and require their continuity at 
the contacts as well as that of u v d x Q v / 'K v . In particu- 
lar, both the charge and spin current, j v = ^/2dt<& v /iT 
are conserved. This is because we neglect interactions 
that violate the conservation of j p and j a , corresponding 
respectively to the umklapp process and the backscatter- 
ing of electrons with opposite spin (g± is irrelevant for 
repulsive interactions, see ref. |l^ for the ITLL). 

In view of the above continuity requirements, the 
charge and spin of the incident electron are reflected at 
the first contact with two different coefficients similar to 
eq.(2.17) with a subscript v = p,cr, 



\-K v 
l + K v 



(3.2) 



The transmitted charge and spin propagate at differ- 
ent velocities u p ^ u a and thus reach the second con- 
tact where they get partially transmitted at different 
times t v = Lju v . Since the transmitted charge and spin 



11 



propagate at the same Fermi velocity in the right non- 
interacting lead, they will stay spatially separated, at a 
distance vp \t p — t a \ (fig||). Due to subsequent internal 
reflections, we end up with a series of spatially separated 
partial charge and spin spikes on the two non interact- 
ing leads. The transmitted spikes correspond to a com- 
plicated superposition of electron-hole excitations. Note 
that in the relevant case of spin-invariant repulsive inter- 
actions, K a — 1, the spin part is not reflected but gets 
directly transmitted to the right lead, while the charge 
undergoes the multiple reflection process. At times very 
long compared to t p (t a ), the series of transmitted charge 
(spin) spikes sum up to unity. Thus the transmission of 
an incident spin polarized flux is perfect. 




FIG. 2. Dynamic transmission of an incident electron with 
spin up. The charge and spin (hatched) are separated even 
in the noninteracting leads. As an example, we consider here 
Up > u a , K a < 1, K p > 1. 



C. "Phase diagram" 

In an infinite wire, interacting electrons with spin have 
a quite involved "phase diagram" ; here we do not intend 
to study it entirely when the wire is connected to leads. 
We give some qualitative results and then discuss thor- 
oughly the tendency towards the formation of a Wigner 
crystal. As for spinless electrons, we expect the dom- 
inant tendency to extend towards the leads in analogy 
with the proximity effect, but the leads can also inter- 
vene in the competition as will be shown later, especially 
at low enough temperature. 

Let us first discuss the case where superconducting or- 
der can take place, as is the case in an infinite wire with 
K p > 1. One has to distinguish singlet SS, corresponding 



to the spin gap case, K* = and triplet SS, K* = 1. If 
we naively take these values for K a inside the wire, then 
we can show that: 

For K* = inside the wire in order to take into ac- 
count the spin gap, we can show that the tendency to- 
wards singlet super-conductivity holds for K p > 1, but 
the proximity effect towards the external leads shows up 
only for K p > 3. Remarkably, for K p 3> 1, an incident 
electron with spin up is reflected back into one hole with 
spin down, and two quasi-particles of opposite spin, of 
charge unity and moving at|-yplpcity u p are transmitted, 
recalling Andreev reflectionBliS 

For K p > 1, and K* = 1, we can show that triplet 
superconductivity develops and extends towards the ex- 
ternal leads, but the Andreev reflection is more subtle to 
interpret .£3 

In the following, we focus on the situation K p < 1, 
K* = 1, as it should be for quantum wires, and determine 
the dominant CDW in the density-density correlation 
function. 



1. Density-density correlations 

Let us compute the density-density correlation func- 
tion. In order to express the density, we can superpose 
eq. (2.2) for the spin up or down, thus with the addi- 
tional index v = p or <r, then express <!>„ as function of 



This yields 



p(x) = -d x 



E 

m p ^0,m CT 



2im r , 



(3.3) 



where the sum runs over integers m p and m a of the same 



parity,H and $ p = $> p — kpx. Normally only m p = ±m CT 
are allowed. But the other harmonics with m p ^ m a are 
generated by renormalization in the induced density_Lhat 
develops as a response to an external perturbation.!!!] 

The fcrmionic correlation functions can be expressed 
through the bosonic correlation function U v (and its 
dual) for v — p, a as in eq.(2.6) with u(x) and K(x) 
indexed by v. Then the density-density correlation func- 
tion is similar to eq. fl2.4| ) if one adds the subscript p to U 
in the first line, and X m is replaced by X m :TO(T obtained 
from eq. (|2.7|) by the following substitution 



2m 2 U -> m 2 p U p + mlU a . 



(3.4) 



Thus 



X.,, 



Xm.X 



where X mv is the analogous of X m without spin with u 
and K indexed by v. The imaginary part of the Fourier 
transform of X m m „ a t low frequencies, Xm p ,m a is de- 
fined as in eq. (Bjj) . It cannot be factored unless the spin 
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and charge velocities are equal. For studying the domi- 
nant tendency, it is sufficient to consider its local values. 
Up to a slowly varying function of X, we have 

{x, X) = <J Xm p (x,x)y/ Xm„ (x,x), (3.5) 

where Xm v { x ,x) is exactly that computed without spin, 
with the additional index v, in particular there are two 
thermal lengths associated with the different charge and 
spin velocity Lt,v — u v /T, two temperatures Tj, „ = 
Uv/L, and two times (in units of tq) 



U„T 



for a spin (y = a) or charge [y = p) excitation to go from 
x to the closest contact. 

The calculation can be carried out for any parameters 
for charge and sp in, o n the leads or inside the wire, but 
let us specify eq.(3.5) for non interacting leads, and for 
K a = 1. In the high temperature limit, using eq.(2.21), 



Xm P , m >,x)=T m o 



m K„+m a -2 



(tanhT^) m p 



1p 



(3.6) 



In particular, at a distance much greater than Lp p from 
the contacts, \m ,m„ (x, x) is the same as that in an infi- 
nite wire 



X n 



Xx,x)~T m ? 



m n K„+m n 



(3.7) 



For m 2 p K p 



>? a < 2, the (m p , m CT ) component of the 
density-density correlation function is enhanced com- 
pared to the noninteracting case. 

In the low temperature limit, using eq.( 2.22j ) 



m „+ m „ 



! Z Lp 



X . p 



(3.8) 



Note that the superconducting correlation functions 
can be obtained in a similar fashion, one has to distin- 
guish the triplet from the singlet superconducting ten- 
dency, and we refer to |(| for more details. 



2. The suppression of the Wigner crystal 

In an infinite wire, the Ak F CDW, corresponding to 
Wigner crystal, dominates the 2k F CDW for very re- 
pulsive interactions, more precisely for K p < 1/3 and 
K a = 1 as one can inspect using eq.(3.7). We show here 
that the leads can induce also a proximity effect in the 
opposite sense, by suppressing the Wigner Crystal at low 
temperature. This tendency is even more important for 
long-range Coulomb interactions to which case one could 
extend qualitatively the results obtained here. 

Our study without spin has shown that the enhance- 
ment of a 2k F CDW persists with leads. We have also 
observed that the 2mk F CDW are suppressed at low tem- 
perature due to the T dependent term T 2 ( m -1 ) in Xm', 



but this effect intervenes only when the backscattering 
potential is studied on resonance, because the 2k F CDW 
dominates 2mk F for m > 2. 

With spin, one has to compare the 2k F to the 4k F 
CDW, thus the correlation functions Xm p ,m a correspond- 
ing respectively to (m p ,m a ) — (1,1) and (2,0). The cor| 



petition is much less trivial than without leads and isEI 
discussed in detail in appendix ^|. Here we summarize 
the results, illustrated by figs j|| and ||. 

We will focus on the K p < 1/3 case, since for K p > 1/3 
the 2kp CDW is dominant at any temperature. Then 
contrary to the infinite wire where the 4fc^ dominates, 
we show here that there is a crossover temperature T c (x) 
from the Ak F CDW to the 2k F CDW as T is lowered be- 
low T c (x). T c (x) decreases monotonically when x varies 
from one contact to the center of the wire where it reaches 
its minimum value, interestingly given by the following 
power of Tl: 



t c (o) = t; 



\0--Kp) 



(3.9) 



One can check that T c (0) < Tl because K p < 1/3. In- 
deed T c (x) saturates at T c (0) for all x such that a— \x\ ~ 
a, i.e. t x T L - 1. Thus for any T < T c (0), the Wigner 




Wigner crystal 
X(x) 



2kp CDW 



" a x 

FIG. 3. The case K p < 1/5: the Wigner crystal is sup- 




a 



pressed below the crossover temperature T c (x), eq.(S6) 



crystal is suppressed all over the wire and the dominant 
tendency is towards the 2k F CDW. 

In order to discuss the case T > T c (0) we shall in- 
vestigate in detail the behavior close to the contacts. 
There is a second particular value other than 1/3 that 
arises: all the correlation functions at ±a are analogous 
to those in the center of the wire with K v re placed by 
K va = 2K V /{1 + K v ) (the analogue of eq. fl2.18D with the 
supplementary index v = p 1 a). When K a = 1, K aa = 1, 
the limiting value K a , P — 1/3 corresponds to K p = 1/5. 
In particular, as long as K p > 1/5 one has K a ^ p > 1/3, 
thus the 2k F CDW dominates the 4k F CDW at the con- 
tacts at any temperature. 
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Wigner crystal 




-a - Xo ^ a 

x 

FIG. 4. The same as fig.|, but with 1/5 < K p < 1/3. 



This yields different expressions for T c {x) depending 



on K„ 



For 1/5 < K p < 1/3, it is given by eq.flCq), and di- 
verges at ±a; reassuringly, we find that T c (x) is con- 
tinuous at the point xq [eq.(C7)] where it reaches 
Tl from below or above (figjj). 



For K p < 1/5, T c (x) is given by Eq.(C6). T c (x) 



is below Tl for any x, and reaches its maximum 
value at a; T c (a) is similar to eq.( |3.9[ ) where K p is 
replaced by K ap . (fig.|) 

If we consider now any K p < 1/3 and a given tem- 
perature in the range [T c (0), T c (a)] (T c (a) is infinite for 
K p > 1/5), we see that the dominant tendency switches 
from the 4fcp to the 2/cf CDW when one goes from the 
center to the contacts. The point x at which this tran- 
sition occurs is a temperature dependent function x c (T) 
that is the inverse of T c (x). For T higher than T c (a) 
(lower than T c (0)), the Ak F (2k F ) CDW dominates all 
over the wire. 

To conclude, the external leads suppress the impor- 
tance of the 4k F CDW at low enough temperature or 
close to the contacts, thus preventing the tendency to- 
wards the formation of a Wigner crystal. 



D. Backscattering by non— random impurities 

The role of impurities when spin is taken into account 
can be treated in a similar fashion as without spin, by 
doubling the indices as above, leading to different pow- 
erlaws on temperature. Nevertheless, in view of our pre- 
vious discussion, the main additional complication arises 
from the competition between the 2k F and Ak F backscat- 
tering: the leads intervene more strongly if K p < 1/3. 

We will first derive briefly the renormalization equa- 
tions at finite temperature, then give results for the con- 
ductance. 



The coupling to impurities is des crib ed in terms of 
/ p(x)V(x), where p is given by eq^SJB). For any po- 
tential, the partition function can be expanded in term s 
of V' in a similar way as for the spinless case |eq.( 2.25 )], 
but with rrii replaced by a couple of integers (m Pt i, m a _i). 
Then Z describes integer charges restricted to two neu- 
tral cylinders, corresponding to the charge (y — p) and 
the spin (y = a) degree of freedom. The radius of each 
cylinder is given by /3/tt, and its height determined by 
the spatial extension of the potential. Only c harg es on 
the same cylinder v interact via U v (x, y, r), eq.( |2.6| ) with 
the index v on any function or parameter. 

Following similar steps as in the spinless case, we can 
i nfer the leading renormalization equation analogo us t o 
( |2.26|) , where U is replaced according to the recipe ( \S-4 ) , 



dV'(x.m p , m a ) 
dl 



1 

l -2 



2 dUp_ 
" dl 



, dUg 

'~dT 



V. (3.10) 



For simplicity we have omitted the arguments on the 
right hand side. In particular, through U v , V now ac- 
quires an additional spatial dependence. 

The bare value of V(x, m p , m a ) is zero if these terms 
are no t present in the initial density representation, 
eq.(|3.3|). In order to illustrate how they are generated, we 
specialize to the case of one barrier. Then the cylinders 
shrink to circles, and the partition function reads 



n J(n-l)ro 



II V ™p*,™*i eX P 



G?Ti 



2J m ul m vj U 1J (n - Tj) 



The sum runs over all n-tuples of integers m P i and m aj 
with total vanishing sum separately, and the bare trans- 
form of V 



V m 



V'(x)e 



2im p <j> (x) 



(3.12) 



acquires a dependence on m CT during renormalization. 
We can show that the interactions are not renormal- 
ized as for the spinless case. The next-leading term to 
eq.( 3.10| ) is obtained by the contraction of one pair of 
charges on each circle, and is equal toE3 



'y(n P ,n<r)V(n' ri a ), 



(3.13) 



up to non universal prefactors, the sum running over n p + 
n' p = m p and n a + n' a = m a . 

For a potential with extension d much less than L and 
Lt, w e can show that the partition function is analogous 
to eq.( [3.1l| ) when a cutoff larger than d is reached. Again, 
at the location of V the interaction parameters might be 
renormalized before reaching this stage. 

If the renormalization of the local effective interac tions , 
as well as next leading terms can be ignored, eq.( 3.10| ) 



(3.11) 
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can be integrated easily. The maximum cutoff is half 
the radius of one cylinder, i.e. (3/2, where we get the 
renormalized parameters!^! 

V'(x;m p ,m a ) = V'(x)xm p ,m„(x,x). 

y m m (x, x) has been stu died previously, and is given by 
eq.(|3.6|) (respectively (|3.8|)) in the high (low) temperature 
limit. 

For K a = 1, K p < 1, and for any x, V is enhanced 
when the temperature is lowered, and is less enhanced 
when one gets closer to the contacts because the leads 
moderate the repulsive interactions. If the potential is 
weak enough or the wire is short enough, we can per- 
form the renormalization at temperatures below Tl . The 
renormalized V'(x, 1, 1) obtained at a given T < Tl is the 
same for all T < Tl, but this is not the case for all the 
other couples (m p ,m a ) ^ (1,1), for which V'(x, m p , m a ) 
goes down to zero at zero temperature. For instance, 
we refer to our previous discussion of the competition 
between the 2kp and Akp backscattering contribution. 

Now we give the perturbative expression for the con- 
ductance. For a general weak pote ntial V(x), the conduc- 
tan ce ca n be expressed as in eq.(2.33). TZ is similar to 



cq.(2.34) with m standing for a couple of integers m 



wit h the same parity. This substitution applies also to 
eq.fl2.35Q, giving 



V'(x)V'(y) 2imp[M*)-My)] v 
(k F u) 2 Xr ' 



The nonlocal function Xm p ,m c r (x, y) is now more difficult 
to compute due to the different velocities of spin and 



charge. But if we assume u p 



it can be inferred 



from that without spin (see appendix [EM) by replacing 



2m 2 K 



m p K p + m\K a 



When the potential has an extension less than L and 
Lt,v (y = P,o~), only the local value, is required as far 
as the dominant term is c oncerncd.Ej Then the leading 
term is similar to eq.( 2.36| ) with m repl aced by the pair 
( m , m g ) , thus replacing V m by eq.( 3.12| ), and using Eqs. 

For K p < 1/3, we refer to figures |3| and |4[ Above (be- 
low) the crossover temperature T c (x), where x is the loca- 
tion of the barrier, the 4/cf (respectively 2k F ) backscat- 
tering gives the do minant contribution, and can be in- 
ferred from Eqs. (3J3 3~8| ) by taking (m p ,m a ) — (1,1) 
(respectively (2,0)) depending on whether T < Tl or 
T>T L . 



E. Extended disorder: Conductance fluctuations 



For a random potential verifying eq.( 2.37 ), the leading 
term of TZ is given explicitly in figure |3[ Again the case 

we have to distinguish 



K p < 1/3 is more complicated 



thr ee t emperature regimes, separated by Tl than T c (0), 
eg.®. 

We can also compute the variance of TZ, using Eqs. 
.14. 2.37 ). Similar steps to that without spin can be 
carried on. While the powerlaw behavior of Var (TZ) is 
different from that without spin, th e rela tive fluctuations 
are the same, i.e. we h ave s till eq.(2.44) in the low tem- 
perature limit, and eq.(2.47) (with a different function / 
that does not affect the order of magnitude) in the high 
temperature limit. 

Let us give the behavior of Var(TZ). In t he low 
temperature limit T < Tl, the inequality ( 2.43| ) holds 
again, s howi ng that Var(TZ) is of the same order as 
TZ 2 , eq.( [2.44 ), given in figj|. In particular, in the case 
1/3 < K p < 1, 

'e 

For K p < 1/3, this powerlaw holds at lower temperature 
T < T c {0), eq.(||), but 

for T c (0) <T<T L . 

In the high temperature limit T > Tl, restricting our- 
selves to K p < 1, we can retain only the contribution of 
the bulk. Then we find 



,{x,y). (3.14) 



Var{1Z) 



U P ToL _ T ni 2 K l ,+mlK c , 



2 — 



(3.15) 



up to a bounded function of both UfLtu for v = p,o~. 
Again, we have ignored the terms due to first and second 
derivatives of \, and kept only the pair (m 



pi 1 



that 

yields the dominant term. This pair is equal to (1, 1) 
(respectively to (2,0)) for K p > 1/3 (respectively K p < 
1 /3) as is the case for t he do minant contribution to (TZ) . 
Remarkably, the ratio ( [2.47 ) is the same as that without 
spin up to the function /, and is -C 1 in any parameter 
range. 



IV. DISCUSSION AND GENERALIZATIONS 

Different aspects of transport in an interacting wire 
with measuring leads have been treated in this paper. In 
this section, we will compare the effect of backscattering 
on the conductance to that obtained without leads, and 
then discuss the results that can be generalized to smooth 
variations of the interactions at the contacts. 

In an infinite wire with repulsive interactions, as 
is thought to be appropriate for a quantum wire, a 
unique local-barrier is sufficient to suppress transport 
completely.EfEl The finite length is intuitively introduced 
as a cutoff, thus preventing the conductance to vanish for 
wea^eiiojugh backscattering potential or short enough 
wireJffl'EH With noninteracting leads and spinlcss elec- 
trons, we recover the results of refs. |37],|5| over all the 
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temperature range only for a barrier whose separation 
from the contacts is of the order of the wire length L, 
and under the additional condition K p > 1/3 for elec- 
trons with spin. 

If we consider extended disorder, the agreement holds 
for K < (3 + \/T7)/2 without spin, and for 1/3 < K p < 
1 + y/2 with spin (see fig [|) . In this parameter range, the 
localization length we infer from the breakdown of our 
perturbative evaluation of 1Z coincides with that found 
in ref. |3| in the weak disorder limit. 

We now discuss the extension of our results to smooth 
variations. As we mentioned before, the discontinuous 
parameter, model can be relevant for edges states in the 
FQHElLSEI, but is not well justified microscopically for 
quantum wires. Thus it is important to know which re- 
sults can be trusted for a more realistic interaction pro- 
file. 

The detailed behavior of the transmission dynamics 
of an incident electron was given previously for abrupt 
jumps.Q'Ej By inspecting the boson Green function, wc 
can maintain the qualitative conclusion; in particular, 
the transmitted charge and spin stay separated in the 
external leads. 

The correlation functions are directly related to this 
process.li3 All the properties of x found in subsection 
[IB1 still hold. To be specific, assume that the leads 



are non interacting, Kp = 1 and that K decreases slowly 
enough from 1 to its b ulk v alue K insid e the wire. Again, 
for K < 1, [cf. Eqs. (|2.1l|) and (§T§)] the tendency to- 
wards the CDW holds inside the wire, but is reduced 
when one goes to the leads, until it reaches its noninter- 
acting value in the leads at a distance u/ (loK) where to is 
the maximum energy at hand, recalling the proximity ef- 
fect. Analogous conclusions hold for the superconducting 
fluctuations when interactions are attractive, the coher- 
ence length being now uK/uj. 

In the presence of a non-random backscattering po- 
tential, the renormalization equations ( 3.10| ) are general. 
The role of backscattering is determined by the density- 
density correlations. Consider for instance a barrier at a 
point x. When x is at a distance of order L of the con- 
tacts, the role of backscatteri ng is the sa me as that in ref. 
|. For K < 1, U (see Eqs. fl2T2| ) and ( |2.36D ) decreases 
when x gets closer to the contacts. The noninteracting 
leads (where a barrier is marginal) tend to reduce the 
importance of the repulsive interactions when one gets 
closer to them. Using the boson Green function proper- 
ties on the external leads (restricting ourselves to the case 
K a = 1 with spin), the temperature dependence of the 
correlation function corresponding to the 2mkp CDW at 
T <T L has a term T m * Kh ~ r . Thus for K L = 1, this 
saturates for m = 1 but goes to zero at zero temperature 
for to 1. In particular, the usual scaling argument in- 
ferring the low temperature behavior T < Tl from that 
at T > Tp is justified only if the 2kp contribution inter- 
venes at all temperatures, which is not the case in our 
model for spinful electrons with K < 1/3, or for any 



parameter range when the resonance in the presence of 
backscattering potential is achieved. 

For a random potential, all the results can be general- 
ized for repulsive interactions. But the inequality ( 2.45 ) 
is more general than that, and holds even if K varies in 
any way inside the wire. 

When spin is taken into account, and for K a = 1, the 
tendency towards a Wigner crystal is expected to be sup- 
pressed at low temperature, more precisely the dominant 
tendency is that towards the 2kp CDW at any T < T c (0) 
[cq.(|3.9|)]. This is because the dependence on tempera- 
ture of the Akp CDW is controlled by Kp — 1, thus is in 
T 4i ^ L_1 = T 3 , vanishing in the zero temperature limit. 

In the presence of Coulomb interactions, the 
ballistic co pd | UC . ta nce is still independent on the 
interactions .BBlia One can make an analogy with the 
case K p < 1/3 to expect the transmitted charge and 
spin to stay separated in the leads, and the tendency to- 
wards Wigner Crystal to be suppressed by the leads at 
low enough temperature. 



V. EXPERIMENTS 

We now comment on the experimental situation. In the 
high-temperature limit, the ballistic value of the conduc- 
tance of one-channel quantum wire 2e 2 /h was observed 
by Tarucha et aZ. Jl3 .as well as in the remarkable clean 
wires of Yacoby et a/O But the low temperature behav- 
ior is more tricky. 

On the one hand, Tarucha et al. fitted the observed 1Z 
with the intuitive law [T 2 + Tl]^ K ^/ 2 of ref. and 
infer the parameter K p = 0.7. We however find a more 
complex dependence of 1Z on the temperature and the lo- 
cation of the barriers. For instance, if the backscattering 
is merely due to the contacts, the exponent 0.7 would cor- 
respond to the local parameter K ap — 2K p / (1+K p ), thus 
to K p ss 0.5. Of course we don't expect the same combi- 
nation of parameters to appear in the case of smooth 
variations of K, but this serves as an illustration; as 
we discussed in section [V, the effective parameter that 
controls the powerlaw increases when one gets closer to 
the wire. Probably the potential configuration depends 
on each wire, and more precise fit and evaluation of K p 
are required. Even more caution is needed if the poten- 
tial is random. While th e con ductance is self-averaging 
at high temperature (eq.(2.47)), its fluctuations become 
more important at lo w temperature, of the order of the 
reduction TZ, eq.(2.44); then one has not to measure sim- 
ply the absolute value of 1Z, but also its fluctuations. 

On the other hand, Yacoby et alS3 fabricated high 
quality quantum wires using a sophisticated technique. 
The most striking result is that the conductance shows 
perfect plateaus as function of the gate voltage. Their 
value is Ge 2 /h per spin mode and channel where G < 1. 
G depends on the temperature and the wire length, but 
is reproducible for different wires made in wells of the 
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same width. This indicates that the reduction in the 
conductance is not due to potential fluctuations, but 
rather to the backscattering of electrons from the two- 
dimensional gas when entering the wire. In the high tem- 
perature limit, the one-channel wire conductance is e 2 /h, 
(G = 1) in accordance with theoretical results. But G 
decreases when one lowers temperature, or increases the 
wire length, in a way similar qualitatively to the law 
rpK-i or j\-k f or an interacting wire with parameter 
K. The main objection of Yacoby et al. against this 
TLL interpretation is that temperature is scaled by the 
Fermi energy which varies with electron density, in con- 
tradiction with the observed plateaus which are density 
independent. Further, the parameter K itself is expected 
to depend on the density. 

The importance of these two effects is however not so 
obvious to assess. First, energy scales different from the 
Fermi energy might play the role of the cutoff. Second 
the way K depends on the density (i.e. the gate voltage) 
is difficult to estimate: a fully microscopic calculation 
of K would be necessary, including the capacitive effects 
of the gates. We have however also to remark that our 
model is not the best suited to describe the complicated 
geometry of the experiments. The transition between the 
two-dimensional gas and the wire is not easy to describe, 
and is different from the usual adiabatic transition. Thus 
we think model closer to the actual experiment needs to 
be developed. E3 



VI. SUMMARY 

In this paper we have studied the properties of an in- 
teracting wire connected to measuring leads. We have 
obtained a number of new results, mainly concerning the 
suppression of the Wigner crystal and its effect on the 
conductance with impurities, and the conductance fluc- 
tuations. We have also completed our study of the trans- 
mission process by taking spin into account. We have 
shown here that the clean wire acts as a separator of the 
charge and spin of an incident electron, which don't re- 
combine in the noninteracting leads. We can hope that 
experimental progress allows this feature to be observ- 
able by appropriate spin and charge detectors. We have 
also given a systematic computation of the nonlocal cor- 
relation functions and their Fourier transform at temper- 
ature T > w. This has allowed to investigate both the 
dominant tendency in and outside the wire and to study 
the role of impurities. 

As first mentioned in ref. 0, we have explained in more 
detail how the proximity effect shows up. For repul- 
sive (attractive) interactions, the charge density (pairing) 
fluctuations are enhanced in the bulk of the wire, and de- 
crease, while staying enhanced compared to the noninter- 
acting case, when one goes beyond the contacts. But the 
novel effect is that the leads influence the competition 
between different tendencies in the wire at low enough 



temperature. This was shown to be the case for very 
repulsive short range interactions K < 1/3. Normally, 
a tendency towards a Wigner crystal shows up, corre- 
sponding to the Akp CDW. When connected to leads, 
such tendency gets suppressed below a space-dependent 
crossover temperature, and is replaced by the 2Uf CDW. 
In particular, this allows for simultaneous 2kp and Akp 
CDW in the wire in a certain temperature range. T c is 
controlled by interactions, thus its measurement could 
provide a way to measure K . The Wigner crystal shows 
up more strongly for Coulomb interactions; we have not 
studied this case, but we expect the same qualitative con- 
clusions to hold. 

We have also studied the role of a backscattering po- 
tential. We have developed a formal framework to study 
its role for any non translationally invariant interactions, 
and explicitly at finite temperature. We have derived 
the renormalization equations for a non-random poten- 
tial V with any extension, but we expect interactions 
to be renormalized by a nonlocal V. If V(x) is weak 
enough for this effect to be neglected, the conductance 
g was expressed in detail in the discontinuous parameter 
case, g is controlled by the interactions, but is influenced 
by the external leads, especially for very repulsive inter- 
actions because the leads suppress the tendency towards 
the Wigner crystal below T c (x) [cf Figs. We have 

computed the conductance fluctuations, whose behavior 
is nonuniversal. But remarkably, the ratio Var(lZ) /1Z 2 
stays of the same order, whether the spin is taken into 
account or not, or whether K > 1/3 or not. It is ~ 1/2 
for T < T L , and - L T /L < 1 for T > T L . We can gen- 
eralize the main results to a smooth interaction profile 
near the contacts. 

One of us (I. S.) would like to thank B. Doucot, D. 
C. Glattli, T. Martin, D. L. Maslov and A. Yacoby for 
interesting discussions. 



APPENDIX A: FORMAL EXPRESSION FOR 
THE CONDUCTANCE 

We derive the analogue of a Dyson equation for the 
nonlocal conductivity. We begin by spinless electrons, 
taking into account spin will the be straightforward. We 
consider a general Hamiltonian 



H tot =H + V($), 



(Al) 



where H is given by eq.(^J), and V($) depends on $, but 
not on its conjugate canonical momentum II = d x Q/ir. 
Let us first give the general equation of motion for $. 
Using [n(a;), 3>(y)] = iS(x — y), any functional T of <f> 
verifies 



.dT_ 



(A2) 



This allows us to write the system 
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ftn = -H = |s x (^*) 



where 



l(x) 



<9V($) 
d$(x) 



(A3) 



(A4) 



Note that uKI has the dimension of a force that inter- 
venes in the time derivative of j = uKU. Eliminating II 
in eq.(A2), the equation of motion for <E> reads 



A** = I, 



where we denote 



Art = 



i r 



-d. 



uKd. 



(A5) 



(A6) 



For an arbitrary nonlocal interaction n ohon tial U(x,y), 
D xt contains also an integral operatorJl2u3 and all the 
remaining analysis can be generalized to that case. 

Note that without V, we have 1=0, and thus recover 
the equation of motion D x t$ = 0. In the following, the 
corresponding correlation functions carry a subscript 0, 
in order to distinguish them from those computed with 
V. 

Next we derive a Dyson equation for the nonlocal con- 
ductivity I, eq.(A4). For any A depending o n $ in a 
local or nonlocal way, we can show, using eq.(A2) and 
eq.(A5), that 



Art {-iB{t) ([$(*, t),A])} = i6{t) ([l(x, t),A\) + S(ty 



OA 
1 d$(x) 
(A7) 

where 9 is the Heaviside function. The Fourier transform 



with respect to time of eq.(A7) yields 



dA 



D x „{{<S>{x);A))„ = {{l{x);A)) 
where we adopt the notation^ 

/'CO 

{{A{x);B{y))) u = -i / e*" ([A(x, t), B(y, 0)]) 
Jo 



(A8) 



(A9) 



We point out that the correlation functions are defined 



here with the total Hamiltonian H to t, eq.( Al). Similarly, 
we can show that 



D yu ((A;9(#))\=((A;I(y)))u 



dA 



Applying these two identities successively to the Green 
function 

G(x,y,u) = <<<&(aO;*(2/)>L, 



D yu ,D xu G(x, y, u) = D yu 5{x - y) + —/(at, y, oj) 



where / denotes 



AtU) 



((l(x);l(y))) u 



dl(x) 



d$(y) 



(A10) 



In order to invert the relation above, and thus to obtain 
Gj(x,y,uj) where x, y are arbitrary, we multiply its two 
members by Go(y,y,Lo), integrate over y, and use 

DyojG {y,y,oj) = 6(y-y). 

Then we multiply both sides by G(x, x, to) and integrate 
over x. We get 



G(x,y,w) = G (x,y,uj) + 



2ioj 



dxdy 



G (x, x, uj)G (y, y, w)f(x, y, u). (All) 

We can also express this equation in terms of the di- 
mensionle ss no nlocal conductivity (measured in units of 
e 2 /h), eq.(|2~3l|): 



v(x,y,u) = <t (x, y,u) - J J dxdy (A12) 

a (x, x, Lo)a(y, y, u)f(x, y, u). 

Th e con ductance measured in units of e 2 /h is given by 
eq-dlH), 



g = lim a(x, y, w) = a(x, y) 



(A13) 



independent of its arguments, which ensures current con- 
servation. It is also real because it is a correlation func- 
tion of two Hermitian fields. Indeed, for two hermitian 
operators A, B, we have 

poo 

({{A;B)) )*=i / ([B\A^]) = ((A;B)) real. 
Jo 

Besides, 



d 

du) 



POO 

(( A '> B ))u,=o = ll \\ A ' B \) dt is Purely imaginary 
Jo 



Thus lim^^o ((3-(x)',Z(y))) is real while its derivative 
at zero frequency is purely imaginary because X is her- 
mitian. If one supposes cr(x, y) finite, one has 



dx / dy 



{{T{x);I{y))\ 



dl{x) 



= 0. 



Then the zero frequency limit of eq.( A12) yields the exact 
conductance measured in units of e z /h 



9 = 9o 



1 -g / dx / dyf{x,y) 



we get 



where 
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f(x,y) = ^Im-^({X(x);X(y))l 



X m (x,y,u) = -i 



is the zero-frequency limit of eq.(AlO) and go is the di- 
mensionless conductance in the absence of V. If one takes 
spin into account, V can depend on the spin field <& CT and 
its conjugate canonical momentum IL a , both commut- 
ing with the field $ p . This is the case for instance for 
backscattering potential, where th e co ntribution of the 
(m pi m a ) harmonics of the density (3.5) to X reads, after 

integration by parts: J V'(x) exp (2im p $ p + 2im (T 3> (T J . 

APPENDIX B: CORRELATION FUNCTIONS 

We will show how to compute the imaginary part of 
the Fourier transform of the nonlocal correlation function 
( p7| ) at frequency lo < T, i.e. x(%,y) [eq.(|2^)]. 

First we clarify the way analytic continuations from 
Matsubara to real frequencies are done when we regular- 
ize by imposing the imaginary time to obey tq < t < 
[3 — tq. Consider any general op erat ors A and B. Using 
the equivalent of the notation ( |A9| ) for the Matsubara 
frequency iu n , we have 



{e" T °(A(x,it + T Q )At(y,0)) 

-er^(A (x,-it + T )A^(y,Q))}. (B4) 

The derivative with respect to to yields, after changing 
t — > — t in the second term above, 



dXr, 



die 



(it + T )X m {x,y,it + T ), (B5) 



where X m (x,y,z) = (A(x, z)A< (y, 0)). According to 
eq.(|B5|), and after the change of variables t — > irTt, 

r ioo+T_ 

z exp \—2m 2 U (x, y, z)] dz, 



dXr, 



did 



1 



-ioo+T 



(B6) 



((A;B)) iUn = 



0-TQ 



TO 

oo 



d{it) {-e^C^o) ( A ( it + t o )B(0)) 

+e io, n (it+/3-r ) ( A ( it + p _ To ) B (0)) | 



where T = ttTtq, and U is given by eq.(2.6). 

Owing to the analytic pro pert ies of U, we can move 
the integration contour in e q.([B6|) as explained now. Let 
us write the integral in eq.( [B6[ ) as 

rioo+T 

zX m (z), 

ioo+T 

where the spatial arguments are implicit. U is analytic on 
the complex band T < Re(z) < 7r/2, as well as zX m (z). 

Let z = /x + iv. We then have T < sin fi < 1. We can 
show that for v — ► oo, 



X m (z) 



-m 2 K L \v\ 



dte~ Unt {e luJ " T ° (A(it ■ 



(B(-it + T O )A(0))} 



Then we can translate the contour of integration by tt/2 — 
T, thus making it along z = it + n/2. On the other hand, 
we can show that X m (it + tt/2) is an even function of t, 
ro)S(O)) (Bl)thus 

+ 00 

tX m (it + tt/2) = 0. 



where we have deformed the integration contour in the 
complex plane, and used 

(A(0 - z)B(0)) = (B(z)A(0)) 

Let us now specialize to the case 



Finally, the integral 
dX„ 



dh) 



) becomes 

200+7T/2 



2ttT 2 



dze -2m*U{z) 



Hr/2 



A = B^ = 



2im& 



where m is normally an integer but can be any real, and 
consider X mi eq.(2.7) 

X m (x,y,r) = {T T A (x, r)A^y, 0)) . (B2) 

Then we have 

(T T A(x,T)AHy,0)) = (T T At(x,T)A(y,0)) . (B3) 

This is because the transformation: $ — > — $ leaves H 
invariant, but A — > A' , Besides, X m is invariant when 
exchanging x and y because of time reversal symmetry . 
Equation (Bl) becomes, after the analytic continuation 

iuj n — ► uj + id, 



This form is much more convenient than (B6) for three 
reasons: the integrand is real, the singularities are 
avoided, and approximations are easier to make. Thus 
the dimensionless function x, eq.(2.S), corresponding to 
m = 1 can be expressed as 

x i x ' v) = 7^2 / x ( x > y> T ) dr - 

J-ioo+ir/2 

This gives eq.(|2.g|): 



(B7) 



x (x,y) = ■\c(x,y)e- 2G ^' x ^ o) e' 2G{y ^^>\ 



T 



(B8) 



where 



C(x,y) 



d( exp [4G(x iy ,i( + tt/2)} 



(B9) 



All the Xm can be obtained by multiplying G by m 2 
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1. Discontinuous parameters 

Let us give the functkm G in the case of discontinuous 
interaction parametersE£l 

oo 

-2G{x, y, z)/K = h(z,V-) + jh (z, v+) + ^ ^ 

r— ±1 p— 1 

[-f 2p h (z, 2pv L + rv-) + -f 2p+r h (z, 2pv L + rv+)] (BIO) 



with 



and 



h(z,v) = — log [sin 2 ttTz + sinh 2 v\ , 



V— = 7rT \x — y\ ju 

v+ = ttT (L — \x + y\) /u 

vl = ttTL/u. 



(Bll) 



(B12) 



7 is given by eq.(2.17). It is worth noting that one can 
deduce G from o~{x,y,uj) as explained in ref. 16, and that 
a(x,y,cj) is related to the transmission process.!] 



a. High temperature limit 

Let us first consider the limit T ^> Tl. We can show in 
detail that a good approximation of the integral C(x, y), 
eq.(B9), can be obtained if we approximate G by 



- 2G(x, y,z)~K {h (z, «_) + 7 [h {z, v+) ~ v+ + log 2]} . 

(B13) 

There is no coherence between the two contacts sepa- 
rated by L ^> Lt, so that only the role of one reflection 
matters. The integral can be expressed as 



C(x,y) 



2 cosh 2K v- 



l '\ ( — , K, 7 if, — + K a \ tanh 2 V—, tanh 2 IL| 



(B14) 

where F\ is the hypergeometric func tion generalized to 
two variables,c3 and K a given by eq.( 2.18| ). 

In the limit of an infinite wire with constant parame- 
ters, we can set K = Kl, i.e. 7 = 0. Then C(x,y) is a 
function of x — y, or of v [eq.(B12)] and becomes 



C K (v) = B(±,K) 2 Fi (k, K,^+K;~ sinh 2 v 



(B15) 



where B is the Beta function and 2-F1 is the hypergeo- 
metric function of one variable. 



At distances far from the contacts compared to Lt, 
thus for v + 3> 1, we have 

C(x,y)siC K (v-). 



In orde r to e xpress x, eq.(B8), we need also to exponenti- 
ate eq.( |BI3| ) for x = y, and z = tq. Thus one gets x( x i V) 
for any x, y. 

a. Local value For \x — y\ <C Lt, we can check 
that C(x,y) ~ C(x,x). Nevertheless, for e -2G(x,x,r ) ^ 

e -2G(y,y,r ) t(J hold; we must naye > L Thug 

strictly speaking the property ( 2.15| ) (now 

Lmin — Lt 

since we are in the high temperature limit) holds only for 
t x 3> l,t y 3> 1. But this is related to the abrupt jump 
of parameters at the contacts, thus w e don 't expect it to 
hold for smooth variations, where eq.( 2.15| ) can be gener- 
alized. Besides bosonization describes large separations 
compared to utq, thus the variations of K as well as the 
behavior for t x <C 1 are not relevant. 

Recall that t x is the time it takes for a "quasiparticle" 
to go from x to t he c losest contact, measured in units of 
the cutoff t [eq.(C7)]: 



a — \x\ 

UT 



Let us write the local value of x- 



(B16) 



K„ — K 



[T 2 + (f - T 2 ) tanh^TrTi,)] 
2F1 Q, K a -K,K a + 1 tanh 2 (nTt x ) 



The Hypergeometric function 2-F1 is unity at the contact 
(v = 0) and stays bounded and nonvanishing all over the 
wire. At x = a, 



X( 



(B17) 



At x = 0, an analogous expression holds with K a — > K . 
Besides, we can show that££l 



: (a,a) < X (x,x) < X (0,0) 



X' 



for K < Kl, while the inverse holds for K > Kl- 

We note finally that the simplified expression ( 2.21 ) 

wa s obta ined by ignoring the hypergeometric function in 

eq.(B16), as well as constant factors, and taking t x ^s> 1. 

But we c an ch eck that x( a , a ), e q-(B17), can be recovered 

from eq.( [2.2l| ) by letting t x = 1. 

b. Expression of x{ a ,x) Another interesting expres- 

sio n is t hat of x{ a , x ) f° r 1 "C t x <C T~ , thus w_ -C 1 in 

cq. (|B14| ): 



X (a, x) ~ \B (\,K a \ Tm«-D tx K a -K (Bl 8) 
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b. Low temperature limit 

The com putation of C(x,y) is now much easier. Since 
in eq.(B9) z = iQ + 7r/2, — sinz = cosh£ > 1 ^ vl — 
smhvL. Then one can approximate, in eq.(BK), 



2G(x, y, < + tt/2) ~ -2K L log [cosht] , 



thus eq. (|B9|) reduces simply to 
dt 



C(x,y) 



(B19) 



(B20) 



independent of x,y. C{x,y) is the same as C(x,x) eval- 
uated in an infinite wire with parameter K^. But x 
[eq. QBSj )] has also parts d epen ding on x and y due to 
G. By letting x — y in eq.( B10| ) we get 



T -l e -2G{x : x,r ) 



(B21) 



up to a bounded and slowly varying function of ^.In- 
serting this expression for each of x, y as well as eq. (|B2C ) 
in eq.(fB"8|), we obtain x( x >y)i factorizing as in eq.(2.1S) 
because L <c Lt- 

Note that taking t a = in eq.(B21) or ignoring the 1 
in front of t x then taking t x = 1 gives the same result 
for x( a i a )- The behavior for t x <C 1 is not very relevant 
since it depends closely on the unrealistic abrupt pro- 
file of u, K, and also because the TLL model is valid at 
distances much larger than utq, in particular t x ^> 1. 



\{x) 



Xi,i 

X2,0 



T}- 3K (tanhTi^ 



K 



(CI) 



where 7 is given by eq.(pT^), and t x = (a — \x\)/u p tq is 
the time it takes for a charge excitation emanating at x 
to reach the closest contact. The Akp CDW dominates 
the 2kp one wherever \(x) <C 1. Thus: 

• For Tt x ^> 1, the Akp CDW is dominant for any 
K < 1/3 as in an infinite TLL, because tanhA ~ 1 
for A > 1. 

• For t x ~ 1, i.e. at the contacts, \(a)T 1 ~ 3Ka : the 
Akp dominate for K a < 1/3, i.e. for K < 1/5. 

• For any t x T <C 1, we can replace tanh e ~ e if e <C 1 , 



and the ratio (CI) becomes 

X = T 1 - 3K H- 3 ^ K . 



(C2) 



This expression holds for t z > 1, but the compu- 
tation can be carried out for any x (see appendix 

Since if < 1/3 < 1, 7>0, t~ 3lK increases when 
one gets closer to the contacts, being <C 1 for t x ^S> 
1, and of order 1 at the contacts. This means that 
the importance of the 2kp compared to the 4fcp is 
enhanced as one gets closer to the leads. We have 
to discuss whether K<l/5orK>l/5. In the 
former case, X(x) <§; 1 (cf eq.(|C2"|)). In the latter 
case A reaches unity for 



t c (T)=T- 



(C3) 



APPENDIX C: THE SUPPRESSION OF THE 
WIGNER CRYSTAL 

Here we discuss in detail the competition between the 
2kp CDW and the 4fcf CDW, separating the high- and 
low-temperature regimes. For simplicity, we restrict our- 
selves to the case K a — 1, and we drop the index p from 
all the parameters or length scales. Note that the effec- 
tive local parameter at the contacts, eq.( 2.18 ), for the 
spin degrees of freedom is also equal to one: the cor- 
relation function U a (x,y,T) is now equal to that in a 
noninteracting wire for any x and y. 



c. High temperature limit 

We consider here T ^> max(TL, p , T^m). Recall that 
the correlation functions at distances far from the con- 
tacts compared to (Lt p , Ltv) are identical to those in 
an infinite TLL, eq.(3__7). At t he co ntacts, the equivalent 
TLL has a parameter K a , eq.( 2.1S ). 

For any x, let us use eqs.( 3.5. 2.22 ) to obtain the ratio 
between the 2fcp and 4kp CDW correlation functions, 



Thus, for t x < t c (T), the 2kp CDW dominates (see 

To summarize the high temperature regime, the Akp 
(2kp) CDW is dominant all over the wire for K < 1/5 
(K > 1/3). In the intermediate parameter region, i.e. 
for 1/5 < K < 1/3, the 4fc F CDW dominates only for 
poi nts far enough from the contacts, i.e. for t x > t c (T), 
eq.©. 



T < 



d. Low temperature limit 

Consider now the low temperature limit, 
inf (7l,p, Tl.o-)- Let us write again eq.(3. 



Xr, 



r (x,x) 



T 



m„*/,m*fK 



L,p 



t x 



(C4) 



This holds for t XtP 3> 1, and the value at the contact can 
be recovered by simple substitution t x = 1 in eq.(Cj 
as we can check from explicit computation for any x 
The essential difference between xi,i & n d X2,0i due to the 
leads, is through the temperature dependent term T 2 in 
%2,o that vanishes in the zero-temperature limit, while 
Xi 1 is temperature-independent. Then the 2kp CDW 
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dominates the 4/cp CDW in the zero-temperature limit 
all over the wire for any parameter K. For K > 1/3, 
the 2kp dominates everywhere at any temperature. For 
K < 1/3, the 4fcp CDW dominates in the high temper- 
ature limit. Thus there is a crossover temperature we 
determine now. 

Let us write the ratio of the corresponding correlation 
functions at finite temperature, using eq.flCMj), 



K\ 



\{x) = — 

X2,0 



L V : K (T L t x ) 
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(C5) 



Again, we have to distinguish the cases K < 1/5 and 
K > 1/5: 

• For K < 1/5, we can easily get the crossover tem- 
perature by letting X(x) ~ 1, 



(C6) 



We can check that for any x G [—a, a], T c < Tl, 
and T c is a decreasing function of t x (fig.^)) . 

For 1/5 < K < 1/3, the same expression for the 
crossover temperature holds, eq.(C6). But it is in- 
side the actual range of temperatures, i.e. T c < 
Tl only for points less than xq where xq verifies 
T c {xo) = T L , i.e. given by 



j. rji 3-yif 



(C7) 



Let us now analyze the high and low temperature results 
together. 

• For K < 1/5, the 4fc F CDW dominates for T > 
T c (x), ec[. (pq) , while the 2k p CDW dominates be- 
low (fig.|3|). In particular, for a given temperature in 
the range [T c (0), T c (a)], there is a spatial crossover 
between the two tendencies occ urring at x(T), one 
can obtain by inverting eq.(C6). 



For 1/5 < K < 1/3, we can invert the relation (C3) 
to get the crossover temperature higher than Tl for 
|x| > Xq, and we can (which is reassuring) check 
that the p oin t where it reaches Tl is again given 
by Xq, eq.(C7), showing that we have a continuous 
function T c (x) given by 



T,lr) = T L [±- 



= TV — 



-h K 



for t x > t X0 (C8) 



for t x < t xo . 



Note that the 2k p CDW dominates for any tem- 
perature at the contacts: t x ~ 1, thus T c (a) in the 
second line diverges (fig|§). 



Contacts, 2kp 



1WJ 




rp-7 




Bulk, 2\ 



T K 7T L 



T(0l^Bulk, 4hp 



r K 7T L 



T L T 

FIG. 5. Tl = 1 - g/(2e 2 /h) multiplied by 1 £ /ut for ex- 
tended disorder as function of the temperature (the x-axis) 
and K p (the subscript p is dropped for simplicity). In each 
of the three regions, we indicate whether the dominant con- 
tribution comes from the bulk or the contacts, and from 
the 2/cf or the 4&f backscattering. The impurities in the 
bulk (controlled by K p ) dominate those closer to the con- 
tacts (controlled by K ap ) either at any T if K p < 1 + y/2 or 
at T > T b = T~ 1/hpKp) if K p > 1 + y/2. For K p < 1/3, 
the 4fc_F backscattering dominates only at high temperature 
T > T c (0), eq.(P). 
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